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A MICROSCOPIC SPIKING NEURONAL NETWORK FOR THE 
AGE-STRUCTURED MODEL 

CRISTOBAL QUININAO 


Work in progress - Preliminary version of March 3, 2015 


Abstract. We introduce a microscopic spiking network consistent with the age-structured/renewal 
equation proposed by Pakdaman, Perthame and Salort. It is a jump process interacting through a 
global activity variable with random delays. We show the well-posedness of the particle system and the 
mean-field equation. Moreover we show the propagation of chaos property and we quantify the rate of 
convergence under the assumption of exponential moments on the initial data. 
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1. Introduction and main results 

In a series of remarkable papers, Pakdaman, Perthame and Salort (PPS) [8j El [10] introduced a ver¬ 
satile model for the large-scale dynamics of neuronal networks. These equations describe the probability 
distribution of the time elapsed since the last spike fired as an age-structured nonlinear PDE. Inspired 
by the dynamics of these macroscopic equations, we propose here a microscopic model describing the 
dynamics of a finite number of neurons, and that provides a realistic neural network model consistent 
with the PPS model, in the sense that in the thermodynamic limit, propagation of chaos and convergence 
to the PPS equation is proved. 

In our neuronal network model, the state of each neuron i is described by a K._|_-valued variable 
corresponding to the time elapsed since last discharge. Of course, this approach is quite different from 
classical literature, where the key variable is the voltage: this is an important originality of the PPS 
model. Neurons interact through the emission and reception of action potentials (or spikes), which are 
fast stereotyped trans-membrane current. In the PPS model, the spiking rate essentially depends on 
the global activity M of the network. Specifically, a neuron with age x (duration since it fired its last 
spike) fires a spike with an instantaneous intensity a(x, M) where M is the global activity of the network. 
Subsequently to the spike emission, two things happen: the age of the spiking neuron is reset to 0, and 
the global variable M increases its value by an extra value of e/N. The coefficient e represents the 
mean connectivity of the network. When no spikes occur, the global variable M is supposed to decay 
exponentially to 0 at a constant rate a. 

For each A € N, let us consider a family (A/”/,..., A/)^)t>o of i.i.d. standard Poisson processes. 
Let us also consider a family (ri,..., tn) of i.i.d. real valued random variables with probability law b. 
These coefficients represent delays in the transmission of information from the cell to whole network. 
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Furthermore, we assume that the family of delays is independent of the Poisson processes. Typical 
examples are 

b=5r{ds), r>0 or b=ce~‘^^ds, c > 0. 

A special case is when r = 0, in that context the network is said to be in a no delay regime. 

Throughout the paper we assume chaotic initial conditions, in the sense that the initial state of the 
neurons are independent and identically distributed random variables. Therefore, for gQ and mo two 
independent probability measures on R+, mo)-chaotic initial states consists in setting i.i.d. initial 
conditions for all neurons with common law equal to go, and setting independently, for the global activity, 
another random variable distributed as mo- 

Our aim is to understand the convergence of the ]R+-valued Markov processes 

fYN\ /’vi.Ar 

Jt>0 — 1^ )t>0, 


solving, for each i = I,..., N and any t > 0: 

ft 


( 1 . 1 ) 




t- / X 

Jo 


i,N 


{du,ds), 


with the coupling given by the global variable 


( 1 . 2 ) 


= M, 


N 




Mf__ )} 


M^[du, ds) 


where is the shifted (in time) process Aft-Tj extended by 0 for negative values of the time. These 
processes are a consistency restriction on the spiking times: when a neuron j sends a signal at a time 
t > 0, the global variable receives it only at instant t Tj. 

We emphasise that it is necessary to impose that the delays are compactly supported, i.e. on an 
interval of the type [—t, 0] for some fixed r. Otherwise the processes would depend on their whole 
past trajectory. Equations are stochastic differential equations on the infinite-dimensional 

space of continuous functions from [—r, 0] to R+, i.e. on the variables Xt = {Xs,s & [t — T,t]) (see e.g. 
Da Prato-Zabczyk [2]). 

Finally, we make the following physically reasonable assumption on the intensity spike function of the 
system: 

{ a(-, •) is a continuous non decreasing function in both variables, 
a(0,-)=0, a(-,0)>0 

a{x,m) CX3, Vm€]R+, 


and impose a second consistency restriction 


(1.4) (V5 > 0)(3a:J > 0) such that a{x,m) < S, Vm G 1R+ 

representing that, independently of the level of the network activity, a neuron cannot spike two times in 
an arbitrary small period of time. To fix ideas, recall the expression of the refractory intensity function 
of Pakdaman-Perthame-Salort 0 ), given by 

{ a{x,m) = 0 for x G (0,x*(m)), 
a{x,m) > 0 for x > x*{m), 

-£^x*{m) < 0, x*{0) = x*_^, x*{oo) = x*_ > 0, 

where x*{m) represents a refractory period, x\ the spontaneous activity due to noise, and x*_ the minimal 
age necessary to spike. 

For the previous setting, we have directly the 


Proposition 1.1. Under hypotheses la and (HID, let N > 1 he fixed and assume that a.s., 

- 71 /TIV ^ 

max An < oo. Mo < oo, 

l<^<N " 

then there exists a unique cddldg adapted strong M.+ -valued solution , Mf^)t>o to (ll.ll) - (ll.2l) . 
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Under suitable conditions (to be preciser later on) we show that the solution {Xj^)t>o behave, for large 
values of N, as N independent copies of the solution to a nonlinear SDE. Let Yq (respectively Mq) be 
a go-distributed random variable (resp. mg) and Aft a standard Poisson process independent of Yq and 
Mq. Then we look for R+-valued cadlag adapted process (1), Mt)t>o solving for any t>0 

nt nOO 

(1-5) Yt=YQ+t— / Yg- / l{u<a(Y'„_,M«_)}-^(dwjds), 

Jo Jo 

and 

(1.6) Mt = Mq — a f Mgds — e 

'-do ..'o Jo 

Remark 1.2. Let us assume for a moment the hypothesis of instantaneous membrane decay, i.e., a is 
going to infinity. In that case equation (dD writes 

rt 


E[a(Fs-u), -Ms_„)] b{dw)ds 


Mt = e f E[a{Xt-ni,Mt-w)]b{dw), 
Jo 


in particular, Mt is a deterministic function of t. Then, if the probability density of Xt can he written as 
f(t,x)dx, the previous relation is reduced to 

nt nOO 


pt nOO 

M{t) = £ / a{x, M{t — w)) f{t — w,x) dxb{dw). 

Jo Jo 

that f{t,x) solves (at le 
a{x,M{t))f{t,x) = 0, 


Coming back to equation (EU, we have that f{t,x) solves (at least in the weak sense) the PDE 

dt dx 

■= f{t,x = f>)= a{x,XI{t)) f{t,x)dx, 
which is exactly the model that motivates our study. 


The nonlinear SDE is clearly well-posed if we, for instance, make a Lipschitz continuity assumption 
on the intensity function. In order to avoid this simplification, here we try to use the approaches of 
Fournier-Locherbach [3 and/or Robert-Touboul [H]. Nevertheless, we are not able to provide a very 
sharp result and we have to restrict ourselves to the case of bounded exponential moments. We introduce 
the quantity 

WUtWr = sup{|Ut|, 0<t<T}, 

which is well defined for locally bounded processes. The second natural result of the manuscript is 

Theorem 1.3. Let us assume that hypotheses (ll.3l) - (ll.4l) hold, then there exists a weak solution (Y), Mt)tyQ 
to (ll.5D - (ll.6D such that 

rt 


(1.7) 


f f E[a{Ys-w,Ms-w)]b{dw)ds 

Jo Jo 


< oo. 


Vt > 0. 


Moreover, if the law of (Yq,Mq) is compactly supported, then there exists a unique strong solution 
(Yt, Mt)t>o to (ll.5l) - (ll.6l) in the class of functions such that there are deterministic locally bounded func¬ 
tions A, B : R+ I—>■ IR+ such that a.s. 

(1.8) \\Yt/Ait)\\T < oo, \\Mt/Bit)\\T < oo, VT > 0. 


The uniqueness result still hold true if the initial datum has a fast decay at infinite. More precisely, 
let us consider the growing restriction 

(1.9) (3 ^ > 2) (3 0 < p < 1) (3 Cj, Cp > 0) : CpX^ < a{x, m) < ^(l 3- x^~‘^ m^“^), 

and suppose that initial condition has bounded exponential moments 

(1.10) < oo, w > 0. 

Theorem 1.4. There exists a unique strong solution {Yt, Mt)t>o to (ll.5D - (ll.6l) in the class of functions 
of locally bounded exponential moments (11.101) . 
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Finally we analyse the chaoticity of the system. To do so, a few more notations must be introduced. 
We denote by ©(R^) the set of cadlag functions on endowed with the topology of the convergence on 
compact time intervals. By definition, each pair {X^’^, belongs to ID)(R^), and then the sequence 

of empirical measures 

N 

2=1 

is well defined and belongs to P(ID)(R^)). 

Note that we intentionally use a two dimensional variable in the dehnition of by repeating 
in each pair which provides one way to deal with the issue of exchangeability^. The third and last main 
result of the manuscript is 

Theorem 1.5. Let us assume that hypotheses (ir3l)-(IIl hold, and that the law of (Yq, Mq) is compactly 
supported, then the sequence of empirical processes p,Nit) converges in distribution to the law of the unique 
process {Yt, Mt)t>o with {go, mo)-chaotic initial states solution to (ll.5l) - (ll.6D . 

If the initial datum has a fast decay (in the sense described in Theorem and if moreover there is 
a positive constant Cq such that 

(1.11) |a(a;, to) — a(a;^ m^)| < Cq [a(a;, to) A a(a;^ to^) |a; — a;^| + |to — to^|] , 

for all x,x',m,m' G R+. Then the convergence of fi]\[{t) remains true. 

In the weak connectivity case, i.e. e G [0,e:o) for Eq small enough, hypothesis (11.111) can be replaced by 

(1.12) |a(a;, to) — a(a; +/i,TO + h)| < Coa{x,m)h, 
for all X, TO G R+ and any h G [0,1]. 


Remark 1.6. 

• Let us notice that condition dnn) is relatively restrictive: indeed, we are asking that the rate 
function a to be Lipschitz on its second variable, and that the fluctuations on the first variable 
are somehow bounded. Indeed, a direct consequence of (11.111) is that 

da{x,m) a{x + h,m) — a{x,m) 


dx h 

The model is conceived having in mind an intensity function of type 

a(x, to) = x^ + l{ 3 ,;^< 3 ;}(aTO -T b) G N, a, 6 G 1 
that formally speaking makes all hypotheses true. 


< Cq a(x, to). 


Mathematical overview. As we already said, the aim of the present work is to give a new microscopic 
point of view of the age structured equation considered in Pakdaman-Perthame-Salort 0111 uni- Therein, 
the model is proposed as a reinterpretation of the well known renewal equation and the microscopic 
derivation is omitted. In Tanabe-Pakdaman m and Vibert-Champagnat-Pakdaman-Pham [TB] authors 
propose a particle system but the question of convergence and chaos propagation is not addressed either. 
Nevertheless, the questions of existence of stationary solutions for the PDF and the numerical/simulation 
aspects of both: the particle system and the limit equation, are deeply studied and several very interesting 
results, regarding the existence of oscillatory solutions, are given. Moreover, the effects of the finite size 
of the populations are contrasted with the solutions of the limit equation. 

The specific mathematical tools used in the present work can be easily traced down to two recent 
manuscripts addressing the question of chaoticity of a unidimensional model: Founier-Ldcherbach [5] and 
Robert-Touboul m- The first paper solves the problem under the merely assumption of integrability on 
the initial condition, which is a remarkable weak hypothesis. Nevertheless, the path-wise uniqueness proof 
(which is at the end the key point of the method) uses a particular distance that is closely related to the 
equation itself. A different approach, based on the discretization of the limit equation, is presented in 0 . 
There, the convergence is proved imposing compactness on the support of the initial conditions. In |12j . 


^One could think of other ways to deal with this issue. For instance, we could define a set of auxiliary variables 
each of them solving the equation II1.2II with different i.i.d. initial conditions. 
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the authors also assume boundedness of the initial conditions, and provide a qualitative characterization 
of the qualitative properties of stationary solutions and their stability in the finite-size and mean-field 
systems, and a detailed discussion of bifurcations, stability and multiple stationary solutions is given. 

The specific mathematical tools used in the present work can be easily traced down to two recent 
manuscripts addressing the question of chaoticity of a unidimensional model: Founier-Locherbach and 
Robert-Touboul |12j . The first paper solves the problem under the merely assumption of integrability on 
the initial condition, which is a remarkable weak hypothesis in comparison with other models. Neverthe¬ 
less, the path-wise uniqueness proof (which is at the end the key point of the method) uses a particular 
distance that is closely related to the equation itself. A different approach, based on the discretisation 
of the limit equation, is presenting in De Masi-Galves-Ldcherbach-Presutti [3] . There, the convergence is 
proved imposing compactness on the support of the initial conditions. 

To prove the chaos propagation property there are two classical approaches: whether we use the 
coupling method or an abstract compactness argument. The coupling is very intuitive and apply in a 
very wide range of applications. Nevertheless, usually it is assumed that functions involved are Lipschitz 
continuous (see e.g. m)- Moreover, the method provides the rate of convergence by explicitly estimate 
the difference between the empirical measures. In Bolley-Cahizo-Carrillo [1], it is proved that the method 
still apply in the case of locally Lispchitz continuity, but imposing some exponential moment conditions. 
The second method is much more general and uses a more abstract framework. It was introduced by 
Sznitman [S], and it is useful to prove also the existence of solutions to the SDE, but do not provide the 
rate of convergence. 

The present model has two main novelties. The first one is that the system is not one-dimensional 
but add an extra equation for the coupling variable. This issue implies in particular that attempts to 
use the distance of Founier-Locherbach [5] fail unless a more suitable distance is found. Moreover, the 
two-dimensional nature of the empirical measures makes that the rate of convergence attained for a 
Wasserstein distance is lower than A perspective of the work is to use a combined PDE/SDE 

approach to find a sharper entropy function that allows us to have the necessary uniqueness result (see 
e.g. Godinho i)- The second novelty of the present work is the presence of delays, that is central in the 
original PPS model. This issue is solved by using the independence of the random environment (see e.g. 
Touboul [llllTT]), moreover, the extra terms are treated as locally square integrable processes getting the 
sharp convergence rate expected. 

Plan of the paper. The paper is organised as follows: Section [5] deals with the well-posedness of the 
particle system by finding some nice a priori bounds of the solutions. Sectionj^is related to the path-wise 
uniqueness question of the mean field system. The result we get remains very restrictive and only can be 
applied for the case of compactly supported initial conditions and/or fast decay at infinite. Well-posedness 
of the mean field system is studied in Section IH in particular, we use the compactness argument to find 
that the sequence of empirical measures converges to a weak solution of the SDE. Finally in Section[51 we 
use the coupling method to find different rates of convergence. An appendix completes the present work 
with some general well known results in stochastic calculus theory that are useful to our developments. 

2. Study of the particle system 

Throughout the present section we fix the number of neurons A^ > 1. For g, € P(D(R^)) and ip G 
C'(R^), let us denote the duality product, 

{p{t),ip) := [ (^(7t,/3t)p(d7,d/3). 

./D(Rd) 

Let us notice that for any fixed random environment, i.e., any realisation of the initial conditions, the 
Poisson processes and the delays at t = 0, we can construct explicitly a unique solution to the problem (it 
suffices to arrange the jumping times and construct the solution by solving the equation in between two 
consecutive jumps). Therefore there is a unique strong maximal solution [X]’^ ,..., , M^) defined 

on a time interval of the type [0, r*), where r* is given by 

T* := inf {t > 0 : min(|A/’^|,..., \X^’^\, |M/^|) = oo}. 
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Proof of Provo sition \l.l[ If we are able to prove that r* = oo, then the particle system (I1.2|) - (I1.1D is 
locally strongly well-posed proving ProDOsition ll.il Then it suffices to find some nice a-priori bounds. 

We start by noticing that any solution ,..., , Ml^)t>o to (ll.ll) - (ll.2l) satisfies a.s. 


( 2 . 1 ) 


max Xf'^ < maxX^'^ 1, Vt > 0. 

l<i<Ar * “ u J - 


Moreover, denoting by X^ the (empirical) mean of {Xl’^), it follows that 


1 ^ 

"'O Jo 


and using that X[^ is nonnegative, 


1 

N 


pt POO 


Next, we fix some ^ > 0, and use the consistency condition (m to get that 
(2.2) (3a;|) (Vm £ ]R+,x < x}), a{x,m) < S. 

In particular, using that a; > (1 — lx<xj), we have 

^ pt poo 1 ^ ^ pt poo 

N 


1 ^ /*t poo ^ 1 ^ 

i—1 i—\ ^ ^ 


and taking ^ < 1, we find a positive constant Ci such that 

^ pt poo 

ri,N\ 


(2.3) 


where 


1 

N 


Pt pOO 

E / (1 + Xl’^) / ds) < C,{t + X^) + , Vt > 0, 

i=l JO JO 


N pt 


_ Pt pOO 

Zf:=iV-iE/ / l{u<i} ^fHdu,ds). 

Jo Jo 


Finally, from equation (O 


rrF ^ i-t-Ti 

< Mo^ + ^E 


j^,J0 JO 


*-{ti<a(XfL^.Mj^)} X'^{du, ds), 


but a is nonnegative, therefore 

< Mo^ + — E / / )} ^du, ds) <M^ + ae (Ci(t + X^) + Zf), 

j—\ *^0 

As a consequence, a.s. the following estimate holds 

(2.4) + C[{t + X^ + Zf), Vt > 0, 

for another positive constant C[, depending only on Ci, a and e. 

Summarizing, inequalities (EID and (1^ implies that a.s. the solutions to (fT^-dni) are locally 
bounded, and therefore r* = oo. □ 

Let us now study the integrability of the intensity rate a. To that aim, let us recall some general results 
of stochastic calculus for jump processes. First, we notice that for any 1 < i < A^ fixed, the particle 
system equation can be rewritten by 

pt poo 

ri,N 


XlX = xli^ + t- I X 


s- I '^{u<a(X'^’’^ ,M’^ )}d^ ^du,ds), 

0 Jo I - I 
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pt nt poo 


^ N ,t 
OL £ V—^ 

j¥^i- 


0 ^0 


^{s<t-rj}^{u<a(Xi'^ ,Ms_)} (^.U, ds), 


here we see clearly that when the process A/”* has a jump, then both variables are changed at the same 
time, unless the jump is in the time interval (t — Tj, t]. Let f € (R2) a regular test function, thanks to 

the previous remark, the Ito’s formula writes 

(2.5) ) = /(X*’^, M^) + f d^fiXl’^, X/f) ds-a f X/f ) ds 

Jo Jo 

+ J^ [/(0,Mi^ + ^l{,<,_,,}) -/(X:’^,Mj^)]a(X:-^,Mf)ds 

If we look carefully, in the no delay case (r^ = 0) equation (I2.5|) is coherent with a instantaneous jump in 
both variables. The respective local martigale (XIj^(s)) is therefore defined by 

(2.6) M^it) := f [/(0,Xff + - /(X:L^,MiI)] 

J 0 


E 


f) - /(^sLE,^i!+D-) 


'-Jo 


- a{Xl’^, X/f) ds 


Lemma 2.1. If the intensity function a(-,-) satisfies the growing condition (dii then, any solution 
{Xl '^,..., X^’^, to (I1.1X(I1.2I) with initial data such that 


(2.7) 

with q € {1,..., 5}, satisfies 

( 2 . 8 ) 


max E[a(Xg M,)^)'^] < oo, 


l<i<A 


E[{^N{t),a‘^)] < oo, t>0. 


In the no delay case, there is Eq such that in the weakly connected regime e G [0,eo), the growing condi¬ 
tion (|l.lll) can be replaced by (11.121) . 

Proof. We start by applying Ito’s formula (12.51) to /(•,•) = a(-, •)'? (using a stoping time if necessary) to 
get 

aiXf^,M,^r = aiX^f^,M^y + q f a{Xl’^, d,a{Xl^^, ) ds 

Jo 

-qa f x/f a(x:’^, Mf )«-i d^a{Xl’^, X/f) ds 
Jo 

+ 1^ [a{0,Mf + - aiXl’^,M^y] a{Xl^^,M^)ds 


+ E / +jpy- ^ ) ds + M^fit). 

j¥=i 
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Using that a(0, •) = 0, and that a is non-negative and non-decreasing function in both variables, we get 
that 

< a(X^’^,Mo^)^ + gCo T Mf)« ds 

Jo 

+ 9 E r ) ds + (t), 


implying that 

a(X:’^, )« < + qCo f aiXl’^, Mf)« ds - f a{Xl’^, Mf )«+i ds 

Jo Jo 

+ qC^ae f a{Xl^^,M^Yds + qCo^Y. f a{Xi^^, M^Y ds + {t), 

Jo ^ ^^Jo 

Next, we multiply by N~^ and add over z € {1,..., iV}, to get 

{qLN{t),aY < {d'N{0),aY + qCo{l + e) f (^Ar(s), a«) ds 

Jo 


taking expectation and thanks to Gronwall’s lemma we get the conclusion. 
In the case of (11.121) . coming back to Ito’s formula (12.5|) . we get 


Y 1 ^ ■ 

+ qCoae W(s),a«)ds-f — 

® i _1 


.af(s), 


YK'^^MYY = a{Xi^^,XlYY + qC^ f a{Xl^^ ,XIYY ds- f a{Xl’^ ,XlYY^^ ds 

Jo Jo 

N .t 

+ dCoJ2 ) ds + W, 

j=l do 


therefore 


{^iN{t),aY = {^JLNYJ),aY + qCo f (/X7v(s), 0 *^) ds - f (/X7v(s), a«+^) ds 

Jo Jo 

.t N 

+ 2qCoae / (^Af(s), ds -I- N~^y^MaY jt), 
do fcl 

taking expectation we see that the result holds as soon as 1 > 2 g Cga e. 


3. Path-wise uniqueness of the mean-field system 

The object of this brief section is to prove under what circumstances the mean field equations (IE3- 
(jl.611 are well posed. We restrict our analysis to a very strong hypothesis that englobes many of the 
applications that can be considered. 

We start by stating some equivalent upper bounds on the solution to the limit equation (ini)-(IE2l) 
for the mean field system. The proof is very similar to the arguments used to get (EID and (El and 
therefore we do not go into full details. 

Lemma 3.1. Any solution {Mt,Yt)t>o to (II.Ill - (11.211 satisfies a.s. 

(3.1) Yt<Yo+t, Vt>0. 

Moreover, there is a positive constant C 2 , depending only on the parameters of the system, such that 

(3.2) [ E[{l+Ys)a{Ys,Ms)]ds < C 2 {E[Yo]+t), Vt > 0. 

Jo 
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As a consequence, there is another positive constant such that a.s. 

(3.3) Mt < Mo + C''(t + E[Fo]), Vt > 0. 

Proof. Inequalities and (ESI) are easily checked by recalling (HI. The last one, is a consequence of 
a change of variables. Indeed, we have that 


f f E[a{Ys_^„,Ms-ui)]b{dw)ds = f b{dw) f E[a{Ys-w,Ms-w)]ds 
Jo Jo Jo Jw 

< [ bidw) [ E[aiYs,Ms)]ds. 

Jo Jo 

Finally, we use again that y > ysi^ — '^{y<ys}) to get 


E[a{Ys,Ms)]ds < (^[1))] + t) + / E[a{ys, Ms)]ds, 


and the conclusion follows. 


□ 


Proposition 3.2. Path-wise uniqueness holds true for the mean field system (ll.5ll - (ll.6l) . in the class of 
processes (Yt, Mt)t>o such that there exist deterministic locally bounded functions A, B : R+ i-A K+ such 
that a.s. (lEHl) holds. 

We remark that thanks to inequalities (13.11) and (13.31) we have that if initial conditions are compactly 
supported, then there are indeed some deterministic locally bounded functions A, B as in the previous 
proposition. 


Proof. For any two solutions (Yt, Mt)t>o and (Yf ,Mf)t>o, driven by the same Poisson measure J\f and 
identical initial conditions (Yo,Mo) = (Yq,Mq), it holds 


E[\Mt-Mi\] 


< 


< 


ae / E[\a(Ys-nj,Ms-nj) - a(Yf_.,^,M'^_,,,)\] 

Jo Jo 

ae [ E[\a(Ys,Ms) - a(Yf,M0\]ds 
Jo 


b{dw)ds 


and 

E[|Ft-r/|]< / E[\Ys+Y:\\a(Y,,M,)-a{Y:,M')\]ds. 

Jo 

We know that a is a regular differentiable function, therefore it is Lipschitz continuous and bounded 
on compacts. Since both coordinates are bounded for some deterministic {A{f),B(t)) locally bounded 
functions, it follows that there exists Ct a constant depending only on an upper bound of A and B and 
the time horizon T > 0 such that 


E[|yt - Yf\ + \Mt - M(|] <Ct[ E[|F, - Yf\ + |M, - M’J] ds, VO < t < T, 

Jo 

and the conclusion follows by applying Gronwall’s lemma. 


□ 


These strong estimates are actually a consequence of the fast that assuming bounded initial conditions 
essentially reduces the system to considering Lipschitz continuous intensity rate. We can reduce these 
assumptions by asking fast decay on the initial conditions. More precisely, since 

Yt<Yo + t<R [Yo + t < R} C {Y < R} => P(yo + t<R)< P(rt < i?), 

then, if initial conditions have an fast decay at infinity as (II.IOL then 

F{Yo>R) = 

for some uj positive and f given by it follows that 

E{Yt >R)< P(To >R-t)<C,, 
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and an equivalent inequality holds true for Mt, 

P(Mt >R)< P(Mq >R- C'^{t + E[Fo])) < 

Proposition 3.3. Under mean field condition (|l.lll) and growth restriction (HU, path-wise uniqueness 
holds true for the mean field system (fT3])-(ra. in the class of processes {Yt^ Mt)t>o such that for any 
t > 0 it holds 


sup 

0<s<t 




< oo. 


Proof. We start by noticing that since exponential moments are bounded, then all polynomial moments 
are bounded. In particular, 

E[a{YuMt)'^] < E[C|(1 + < oo, 

for any time f > 0. Thanks to hypothesis (ll.llL we also have that 

|a(n, Ms) - a{Yf, M')| < Co {a{Ys,Ms) + a{Yfi M's)) |W - + Co\Ms - M'|. 

Consider now 

£{Y,M) = { sup Ys < R, sup Ms < i?}, 

0<s<i 0<s<i 

then 

E[|yt-P/l] < 4Coi?a(i?,i?) / E[\Ys-Yf\]ds + 2CoR [ E[|M,-M'|] ds 

Jo Jo 

+ [ E[(W + F;)4]i/4E[|a(y^,M«)-a(y;,M')|4]^/^]P(f-(y,M))i/4p(f“(y',M'))^/^ds, 

Jo 

and 

E[|Mt-M;|] < 2aeC'oa(i?,i?) / E[|W, - y;,|] ds + oeCo / E[|Ms.-M',|] ds 

Jo Jo 

+ ae [ E[\a{Ys>,Ms')-a{Yf,M's,)fY^^F{£‘'{Y,M))^/^F{£''{Y',M'))^^'^ds. 

Jo 

Using Gronwall’s lemma, we conclude that there exists a constant Ct depending only on the parameters 
of the system, such that for R large enough 

E[|yt - Yf \ + \Mt - M;|] < CTe'^^^“(^’'”)p(£:=(y,M))l/4p(£■^(y',M'))^'^^ 
but using the fast decay at infinite 

E[|W - Yf \ + \Ms - M's\] < 

where the constant depends on the time horizon T, but not on R. Finally thanks to hypothesis ()1.9|) . we 
get the conclusion: 

exp (CrRaiR^R) - 0. 

□ 


4. The mean field system & the compactness method 

So far, we know that under relatively weak assumptions on a(-, •), for each iV > 1, and {go, mo)-chaotic 
initial states, there exists a unique solution to (HU-dllll). Now we study the convergence of the this set 
of solutions as N goes to inhnity, i.e. the tightness of the sequence of empirical measures p,N- To that 
aim, we start by recalling (see e.g. Jacob-Shiryaev [3 Theorem 4.5, page 356]): 

Aldous tightness criterion: the sequence of adapted processes {X]'^,Mf^) is tight if 
(1) for all T > 0, all e > 0, it holds 

lim limsup sup P(|M_^ — Mg,\ > e) = 0; 

5^0+ AT-i-oo (S,S')gAs,t 

where As^t is the set of stopping times {S, S') such that 0<S<S'<S-\-6<T a.s., and 
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(2) for all T>0, 


lim sup Pf sup (Mf^ +X}’^) > k] =0. 


= iV>l ^tg[o,T] 


Proposition 4.1. Under hypothesis (O and (HI). Consider two probability distributions gojiu-o such 
that 


/ + rn^ + a{x,rn)‘^)gQ{dx)rnQ{dx) < (X, 


and the correspondent family ,... of solutions to (ll.ll) - (ll.2l) . starting with some 

i.i.d. random variables X^^ with common law go, and independent of the mQ-distributed Mq . Then 


(i) the sequence of processes [X]’^is tight in D(R^); 

(ii) the sequence of empirical measures giy is tight in P(D(R^)). 


Let us remark that the sequence := {Xl'^,M^) is exchangeable, then (ii) follows from (i) thanks 
to Sznitman [B Proposition 2.2-(ii)]). 

Proof. We only need to show the Aldous tightness criterion, let us notice that the second condition is 
easy to show. Indeed, from estimate m we get that 

E[ sup Xl’^] <E[Ao^’^]+T < ( f x^ goidx)y^\T < oo, 

0<s<T ^Jr+ ' 

and recalling (1^ . we notice that 

E[ sup < E[M^]+C[{T + E[X^]+E[Z^]) < oo, 

ie[0.T] 

because Zf^ is the mean of N i.i.d Poisson(T)-distributed random variable. We deduce that the expec¬ 
tation of the lefthand side is finite, independently of the value of N. 

Under the assumptions on the initial distributions, the first point is not very difficult to prove either. 
Indeed, by definition 

pS' p<x> 


and 


\Ms - M^,\ < ..Mf_^.) (du, ds) + a / Mf ds. 


Moreover, using Markov’s inequality 

nS' /•oo 


S JO 


^s- '^{u<a(xl'f’ {du,ds) > 0^ 


< E 


rS+5 




< (5^/^ ^ 


1 / 2 -, 


which is finite independently of the value of N thanks to Lemma 12.11 The first term in the second 
inequality is handled in the same way. Finally, 

rS' . ^ , rT , 1/2 


E 




Xl^ds < E {M^fds')^ 


which is also finite thanks to the assumption on the initial condition. 


□ 


The natural next step, in the proof of existence of solutions to the nonlinear SDE, is to prove that any 
limit point of the tight sequence /r at is a solution of the mean field limit system, which is usually called 
consistency of the particle system. This result is stated in the following 
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Proposition 4.2. Under the same hypotheses of Proposition [7^ any limit point p, of pN ci.s. belongs to 
S-= {ciiYt,Mt) ^>o) • ^ SoluttOTl to (fr5])-(fra such that 

C{Yo) = go, C{Mq) = mo and f f ^[alYs-ui, b{dw)] ds < oo, Vt>o|. 

Jo Jo ^ 

A previous step that simplifies the proof of this result is the 
Lemma 4.3. Let us consider t > 0 fixed and define ttj : D(IR^) —>■ by 

Then, Q € P(D(R^)) belongs to S if the following conditions are satisfied: 

(a) Qottq^ = {go,mo); 

(b) for all t > 0, 

/ / / a{'ys-w,fJs-w)b{dw)dsQ{d'^,dl3) < oo; 

Jd(r^) Jo Jo 

(c) for any 0 < si < ... < Sk < s < t, any ipi,... ,ipk £ C' 6 (M^), and any (p G C'^(IR^), it holds 

F{Q) ■= [ [ Qid'l\dP^)Q{d'y^,dl3^)pi{-fl^,PlJ...(pkill,,l3lJ 

J]D)(R^) ./D(Rp 

- [ dppijl,, z?],) [ - 0 / 3 ], +ae [ a{-;l_^,l3l_^) b{dw)] ds' 

- f , Pi')ds' - f 0 ( 7 ],, /3l,) [(/?( 0 , /3l,) - p{jl,, Pi ,)] ds' = 0 . 

J S J S 

Proof. Let us consider a process {Yt,Mt)t>o of law Q which satisfies (a), (b) and (c). From (a) and the 
independency of mo and go, we have 

•^(^0) = go, F{Mo) = mo, 

and from (b) we have that 


f f 'E[a{Ys-u],Ms-w)]b{dw)ds < 00 , Vt>0. 

Jo Jo 


Finally, from (c) we have that for any (p G (R^), the process 

p{Yt,Mt)-p{Yo,Mo)- [ dpp{Ys,Ms) -aMg + aef E[a{Ys^s, , Ms^s,)]b{ds') 


ds 


-f d^p{Yg,Mg)ds- j a{Ys,Mg){p{0,Ms) - p{Yg,Mg))ds, 
Jo Jo 


is a local martingale. The conclusion follows as an application of Jacob-Shiryaev [3 Theorem 11.2.42 
page 86 ] and [3 Theorem III.2.26 page 157]. This result is classic, but for completeness of the present 
manuscript, we provide some remarks on Appendix]^ □ 

We finish this section by giving the proof of Proposition 14.21 

Proof of Proposition^^ At this point, the problem is reduced to prove that conditions (a), (b) and (c) 
of Lemma 14.31 hold. Since we do not have much information about p we cannot work directly with it. 
On the other hand, we know that pN (up to subsequence) is converging to p, therefore, it seems natural 
to use equations (ll.ll) - (ll. 2 l) adequately and then pass to the limit. 

Step 1. Let us recall that for any A^ > 1, the random variables are i.i.d. with common law go, and 
that Mq is random variable with law mo, independent of {X^^). It follows that 


= ( 50 , mo). 


N 


N—¥oo N 






ON A POISSON COUPLING FOR AN AGE-STRUCTURED NEURONAL NETWORK 


13 


We also have, by the Fatou’s lemma and inequality (12.31) . that 


E 


) Jo Jo 


[a{'^s-w,lis-w) A K\ b{dw) ds fJ.{d'y, d/3) 

^ Mf_J] b{dw) ds < oo, 

A _1 'J 0 'J 0 


for any t > 0. Letting K ^ oo we get (b). 

Step 2. It only remains to prove (c), to that aim, we start by noticing that writes 


N 


*- J s J s 

- f [^(0, Mi7) - Mf, Mjy)] ) ds' 

J S 


-I a{Xi;^^,XlM)Kdw)ds' 


N 


N 


i=i' 


At the same time, using the Ito’s formula (12.51) to the test function (p{-, •), we have 

f d,Mf,M^)ds'-a f M^dpMf.M^)ds' 
Jo Jo 




as 


\ 

j — (p{Xj ,Mj,,)] J {du,ds ) 


pt — Tj nOO 

+ E/ + / l{u<a{Xi^,^,M^)}M^{du,ds'), 

j^ido Jo 

implying, that F/p.^) can be rewritten by 

PM = ^ E Mlf, ■ MXlf, + (Aj’^ - , 


2=1 


with 




pOO 

ds) - a{Xs- , Pd^-) ds 


and 


a: 


i,N _ 


/* ^— r / c" \ 

j^i 0 0 


/o ^0 

rt 


pt ^ ^ ps 

- ^ / a{XiX^,MMKdw)ds. 

Jo -''' Jo 


Using that the Poisson processes A/"* are i.i.d., we get that the compensated martingales R\X are orthog¬ 
onal, and thanks to the exchangeability, we get that 

E[|i^(M^)|] < ^E[(i?i’^-i?A^)2]V2 + C;^E[|Ai’^| + |A^^|] 
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for some positive Cp depending on the upper bounds of the test functions composing F. Moreover, the 
first expectation is bounded uniformly on N: 

rt 


<Cf [ E[aiXl’^,M^)]ds, 
Jo 


which is finite thanks to (12.31) . 

For the second expectation, we split Aj’ in four quantities that can be handled separately: 


I r^~'^3 (=• 

+1E / + ^) - )] 

3 = 1 


N pt — 7 

E 


i=i 


/o 


X 

“-^0 

^l,iV ..N , 




ds 


N 


3 = 1 


+ C'f^I E ( / ) rfs - / / a{Xs-n,^M^-w) Kdw) ds^ 

■i 0 0 


: — 


The first three terms are controllable simply using that Lp G (7^(18.^). Indeed, for any (x, m) G M+, we 
have that 

(/j(x,m+^) = (p(x,m) + ^dpip(x,m) + d^^ip(x,m) + 0(N-^), 

using again (j2.3ll we get that the respective expectations are going to 0 when N goes to infinity (using 
Holder’s inequality to find the convergence). 

The contribution of must be handled more carefully, nevertheless, we have that 


E, 


delays 


a{Xi'ZpM:L,^)ds 


OO pt pt ps 


/ / a{Xl’_^,Mf_Jdsb{dw) = 

h Jw Jo Jo 


a(Xfy^, b{dw) ds. 


for any 1 < j < iV, therefore, each term defined by 

:= [\{Xi’Z,,M^_,^)ds- f [\{Xi’^^,M^_Jb{dw)ds, 
Jtj Jo jo 

has zero expectation. Using that the delays Tj are i.i.d, we get that 

E[.^-, 4 (t) X 5"fc,4(t)] = E[Edelays[^-.4(s)] X Edelays[■5fe,4(s)]] =0, 

if i 7 ^: fc, then 

E[^4] < Cp^(^ f2 X = Cp 

i,fc=i ' 

Finally, using Lemma 12.11 we see that 

E[^i.4(s)"] < 2e[( f a{XlL^,^,Mf_,Jdsy + 

J Ti 

the first quantity is bounded because 

j\{XlLZpM^_Jds < X (^l\{Xl’^,Mf)^dsy^\ 


\/F 

a{XlLl,M:Ljbidw)ds 


t ps 2 

l.iV j.,rN 


0 Jo 


and the second because of 


r r a{XlLl,M^_Jb{dw)ds < f a{Xl^^,XI^)ds, 

Jo Jo Jo 
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summarizing, ^ is also going to 0 with N. 


Step 3. Before passing to the limit we still need to be sure that no mass is added in the discontinuity 
points of the paths, i.e., we need to check that for any t > 0, a.s., 

m({( 7,/3) ■■ A(7,/3)(t) 7^ 0}) = 0. 

The proof is exactly as in 0 Theorem 5-(iii)-Part 2] but for completeness we give some remarks. In 
order to get a contradiction, we assume that there are some b, d > 0 such that 

P[£']>0, with := {/i({(7 ,/3) : max(|A 7 (t)|, |A/3(t)|) > 6}) > d}. 

Therefore, for any e > 0, it holds 

£; C > d}, sup max(|A 7 i(s)|, |A 72 (s)|) > 6}. 

—e,i+e) 


Moreover, Bl is an open subset of ]D>(R^), then 

Kd ■■= {Q e P((R+) : Q{Bt) >d}c P(D(R^) 

is also an open set. Thanks of Portmanteau theorem we get that for any e > 0, 

liminf e Vld) > e V^d) > > 0. 

N^oo 

On the other hand, for N large enough, the jumps in equation (EH) are smaller than b and then the 
problem is reduced to control the size of the jumps in equation (HU, and in particular to show that 


1 ^ 

1P{HN e Bid) < X! /*+' 1 


yN'^{du^ds)>l^ J 


0 , 


^<a(X7y 

which can be easily done using the same arguments of the proof of Proposition 14.11 
Step 4- Now we see that F is a continuous function at any point Q € P(ID)(]R^) such that 

Q({(7,/3) : A(7,/3)(si) = ... = A(j,/3)(sk) = A(j,^)(s) = A('y,/3)(t) = 0}) = 1, 

and 


/ / / a(js-w,/^s-iu}b(dw}dsQ(dj,d/ 3 ) < oo. 

dD(R 5 _) Jo Jo 


.. 7 ) Jo Jo 

Thanks to Step 2 and 3 we know that our limit belongs to this subset of P(D(K^), and therefore 

E[|F(^)|] < lim limsupE[|F(/rAr)| A df] =0. 

if-J-oo Ar_>oo 


□ 


So far we have built a weak solution (Yt, Ms)t>o to the mean field system such that 

/g/p* E[a(Ts_u;, ^(dw) ds is finite for all times. Thanks to the path-wise uniqueness result of 

Section |3] we can go a little further by providing the 

Proof of Theorem \1.SX Proposition 14.21 gave us already the existence of a weak solution such that (HI) 
holds true. Furthermore, if the initial laws {go, mo) are compactly supported, then thanks to inequali¬ 
ties and dSSl), we get that a.s. 

Yt <Yo+t < Sgg+t=:A{t), Mt < Mq-I--I-E jYo]) < Smo (J2{t+ W\Yo\) =: B{f), 

with (respectively Sma) upper bound of the support of go (resp. toq). Proposition 13.21 implies 
that there is a unique solution such that the previous conditions hold true, and this solutions is exactly 
a process with law g, (it suffices to pass to the limit in the associated particle system). □ 

The existence of a strong solution in the case of exponential decay of initial conditions is similar and 
therefore omitted. 
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5. Mean-Field convergence by the coupling method 


In this final section we use the ideas of coupling to prove a quantified version of the convergence of the 
empirical laws towards the law of the unique process that solves (IE3 and (fra . We start by noticing 
that, thanks to Theorem ll.31 there exists a family of stochastic processes 

(5.1) 

such that 






i>05 


nt /•OO 

+ t- 


and 


=M^ -a 


'^{u<a{Y^:f (dn, ds), 

^ V ^ V 

ds-s f r E[a{Y:Ll,M:’^J] b{dw)di 
Jo Jo 


'-JO Jo Jo 

where the initial conditions and the Poisson processes are exactly as described in the introduction. In 
the following, we use the notation rjiq(t) for the empirical mean associated to the exchangeable family 

By the definition of , it follows that 


^ 2 ,/^ _ y-i,N _ 


0 JO 




N^{du, ds), 


then 


(5.2) < - / E 

^0 


\xi'^ - r;-^| a{xi'^, ) A a(r;’^, mi^) 


ds 


E 


{XlX + Y:X)\aiXlX, ) - aiY:X, Ml^) 


ds. 


Since initial distribution are compactly supported, we notice that Y^’^ < A(t) and < B(t) for 

some locally bounded functions A, B independent of N. Therefore, a is bounded and Lipschitz continuous 
in both variables. Then, for any i = 1,..., iV, we get that 

E[|x:-^ - < Ct fm\xl',^ - y:;^\ + \m^ - ds, 

Jo 

for some positive constant Ct independent of N. Similarly, by the definition of , we have 

, Af N t 

(5.3) -Y,E[\M^ - / E[|Mf-M^’^ljds 

2=1 2=1 

N .t 

+ ^ E / E[|a(X:-^,Mf) - a{Y:X,MlX)\] ds 

i=i -^0 

'[| E ( MIX) ds - l\[a{Y:X, XIl'Z.)] b{dw) ds) 


as, 

+ iv’ 

we finally see that using the arguments of Step 2 of the proof of Proposition |4^ we notice that there is 
another positive constant, that we also call Ct, such that 

2 = 1 2=1 ^ 

getting the 

Proof of Theorem 1 1. 51 - (compactly supported case). Gathering (15.2|) and (15.3|) and using the Gronwall’s 
lemma, we get that 

lEE[(|x:’^-r;’^| + |Mf-M:’^|)] < 
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To finish, we apply Fournier-Guillin [H Theorem 1] with d = 2, p = 1, q = 2 + e, to find that there exists 
a positive constant C independent of N such that 

but, since the initial laws (pojTno) are compactly supported, all polynomial moments of the solution are 
upper bounded by a constant independent of N. Using triangular inequality we get that 

log(l + N) 


< Ct ■ 


m/2 


since S has only one element we conclude that locally in time 

as log(l + iV)/ViV. 


□ 


Let us now explain why we cannot conclude using the same technique in the fast decay case, i.e., to 
mimic the path-wise uniqueness proof. Recalling Ito’s formula (I2.5I1 . we get 

E[|X*-^-y;’^|] < rE[|a(X:>^,Mf)-a(y;’^,M:’^)| X \Xl^^ + Y:’^\] ds. 

Jo 

Moreover, thanks to the fast decay at infinite it suffices to define the event 


£i,N — 


( sup \Xl 

^ 0<f<T 


sup \y:’^\<R 

0<t<T 0<t<T 


}> 


to get that 

E[\Xi'^ -Y^’^l] <2CoRa{R,R) [ E[|X*’^ - y;’^|] ds + 2CoRf E[\M^ - ds 


+ CTm,N)"^^- 


Finally, multiplying by ^ and adding on i, implies that 
N rt . N 


(5.4) <2C'oi?a(i?,R) / 1 ^E[|X:'^ - y;’^|] ds 

i=l i=l 

<.t 1 Af 

+ 2CoR - ^E[|Mf - Mm\] ds + Ct max^P(f,V)i/^. 
On the other hand, the difference between and , is controlled by noticing that 
= -a [ {Ml’^ - ) ds 

Jo 

P^ rrv ^ r ^~'^3 

+ Ha{Y^/Z,Ms-w)]Kdw)ds - J a{Y^’^, ds 


SO. X—' 

+ ivS 


N i-t-Ti ^ 


i=i ■ 


i=i ■ 

^ pt — Ti poo 


ds 


so V 

~ iv 2-^ 




/o ^0 


^{u<aiXm,M^)} [■^^idu,ds)- duds]. 
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The first term on the righthand side has a nice structure, and the third one is controlled by 

)|] ds 
Jo 

+2CoRa{R,R) [ E[|Xf’^ - |] ds, 

Jo 

which is equivalent to (15.41) 

The other two quantities are a little bit more delicate to handle, but using exchangeability and recalling 
the definition of ^ 1 , 4 ( 3 ) (see Proposition I4.2L gives 


E 


N 

EQ. 

iV ^ 
i=i 


E[a{Yjll,Mtl)] h{dw)ds - / Ml’^) ds 


0 ^0 


ea 

m/2 




, 1/2 


and for the last contribution, we simply recall that the Poisson processes are independent, to get 


ea 

iV 


E 


AI pt — Ti poo 


J = 1 -^0 


WHdu,ds) - duds] 


< 


ea 

7VV2 


E 


a(Xb^,Mf)2 


1/2 


ds. 


We gather all the previous inequalities to find 


1 ^E[|w:'^ - + |M:>^ - Mf I] 


N 


N 


< Ct max^P(£:'=,^)'/2 + + CRa{R,R) [ ^ ^E[|X:;^ - y;;^| + - M^\] ds, 

do 

but, even if we use Gronwall’s lemma and condition dLlI, we still can pass to the limit i? —>■ 00 because 
of the presence of the extra term on the righthand side. 


6. Further asymptotic analysis 


So far, we have seen that if initial condition are strongly bounded (in the sense of a.s. lie inside a 
compact) then we get a nice rate of convergence, which is somehow sharp, on the norm. 

We explore now, under what circumstances these results still hold true in a larger space. To that aim, 
we use the a priori bounds and the integrability of the intensity rate. Let us start by a result that tell us 
that the intensity function remains with high probability bounded. 

Proposition 6.1. Under the assumptions dni) and dEl, for any positive horizon of time T fixed, there 
is some positive constant Cmf such that 

(6.5) P( sup ma^{{nN(t),a^),{r]N(t),a'^)) > Cmf) 0, 

as 1/N. 

Proof. Let us recall that from Ito’s formula (12.51) we get, for a generic function /(•, •), that the compensated 
local martingale writes 


Af)^(t) = 




0 ■- 


^Jo 








'0 




■-^0 


ds) - a{X^s’^, Mf) ds 
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Using that the local martingale (t) is locally integrable and that a(0, •) = 0, the associated brackets 
processes, in the case f = a^, are given by 


= I a{Xy\M^fds 


' 


0 


^[xY_ 


N 

+D- 


M, 


N 


-+'^3 


a£s2 

TV'' 


a{xi: 


N 

+ '^3 


Ml+rf lis<t-r,}a{X^s'^, ) ds, 


and, for i ^ j, 


1 }a(X:’^,Mf)ds 




E / k^. 




/o '- 


s_+rfc’ -^s.+Tfe 






Now, define = N ^ J2^=i (t)^ it follows that 


E[( sup =E[(Af^.)(t)] = ^E 

0<s<t iv 


N 


)(t) + 2 ^(Aflf 


i<j 


< 


”/V ’ 


where the last inequality is obtained using the same arguments of ProDOsition l2.ll In particular, we get 
that 

P[ sup M^At) > 1] < E[( sup MAt)?] < 

0<i<T 0<t<T -^V 

We come back to (|2.5p in the case / = to get 

{^iN{t),a^) < (/rA(0),a^) + 2(70 / (aGv(s), a^) ds + 7W^2(t) 

Jo 

-|-2(7oae / sup {fiN{s'),a){fj,N{s),a) ds, 

Jo 0<s'<s 


therefore, conditioning on the event supp^^j. < 1 we get 

sup (/IAT(s), a^) < 1 + (7oQ;e + (/lAr(O), a^) + 2 (7o(l + 2 Qfc) / sup {^]s[{s'),a‘^) ds. 

0<s<t Jo 0<s'<s 


Finally, we recall that 

1 ^ 

W( 0 ),a 2 ) = 

and by the law of large numbers we realise that there is a constant CV such that 

P[W( 0 ),a 2 ) >l+EA{xA,Mm < 

and we deduce that indeed (p,Ar,a^) remains with high probability bounded. 












20 


CRISTOBAL QUININAO 


Finally, let us recall that 




rt,N JlyrN\2 


i,N 71 /fi.N \2 


Co f a[Y:’^,Mi^^) j\[a{Y:’Z.M:’^J] b{dw)ds +M^f{t), 
Jo Jo 


and /p E[a(Fg*l^,, b{dw) da < Ct, with an upper bound independent of N. Moreover, 


= faiY:’^,Ml’^)^ds, 
Jo 


which in expectation is bounded thanks to the fast decay at infinite of the initial conditions. Then, from 
the independency of the Poisson processes, we also get 


[ sup M^ 2 it) > 1] < E[( sup M^ 2 {t)f] = ^E 


1 


■ 0 <t<T 


0<i<T 


iV 2 


N 




< 


Or 

N ■ 


□ 


Finally, we have all conditions to quantify the convergence of the empirical measures in the space of 
functions with fast decay at infinity. 

Proof of Theorem I f. 51 - (fast decay at infinite case). Starting once again with the Ito’s formula (12.51) . we 
have that 


^0 *- 

Moreover, defining the event 

£n = \ sup (/7Ar(t),a^) < Cm/, sup (77Ar(t),a^) < Cm/K 

^ 0<t<T n<-KT j 


ds. 


0 < t<T 


we get that 


E 




ds 


From hypothesis (11.91) and Jensen’s inequality, we have 


1 ^ 
ivEix 


i,7V|2(l+p) 




„l-p 


< 


' - N 

i—1 t—l 

therefore, multiplying by N~^ and adding on i, we conclude 


Y,a(Xp.Ml<fP-P < [WW.a*)]'-'', 


Ive 

N 


N 


N 


\XfN_Y^,N\ 




^ 0 

therefore there is a positive constant Ct, not depending on N, such that 

^ < CTP(f^)'/" + 


ds. 


N 
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On the other hand, the difference between and , is controlled by noticing that 


= -a {Ml’^ - Mf) ds 

Jo 


N ^t-Ti 


+ E[a(y;i^, bidw)ds - ^ E ^ ds 


SO. \—^ 

+ 


^ /*i—T,- r 


i=i ■ 


i=i ■ 

^ pt—Tn pCO 


ds 


CCt V / 


W^idu,ds)- duds]. 


The first term on the righthand side has a nice structure, and the third one is controlled by 

rE[|a(y/’^,Mi’^) - )|] ds 

Jo 

< C'tP(£’w)^^^ + Co [ |] ds 

Jo 


+ Co 


E 


(a(Xf’^, x {X^’^ + Is, 


ds. 


The other two quantities are a little bit more delicate to handle, but using exchangeability and recalling 
the definition of £^i^ 4 {s) (see Proposition I4.2L gives 


E 


ea 

'n 


N 

E 

i=i 


E[aiYjLZ,Mifj]bidw)ds- / , Ml^^) ds 


0 ^0 


ea 

'nTJ2 


E[^i4(s) 


21 1/2 


and for the last contribution, we simply recall that the Poisson processes are independent, to get 


ea 

'n 


E 


N ^t — Ti /•OO 


u 


^{u<a(Xi'^,M^)} {Jj^i.du,ds) - duds) 


< 


ea 

iVl/2 


E 




ll/2 


ds. 


Gathering all the previous inequalities leads to 

E[|x:’^ - F;’^| + |M:-^ - I] < ^ + E + + Ct - Mjy|)] ds, 


from where we finally get that 

E[(|X*>^ - + \Ml’^ - p)] < max ( 

we notice that the best convergence rate is in the case of p = 1 / 2 . 


— _L') 

/Vi/2’ iVP/’ 


□ 


Appendix A. General Theorems for Stochastic Processes 


Remarks on the proof of Lemma 14.31 Let us recall that the last thing we got in the main text was 
that for any ip G the process 


(A.l) ip{Yt,Mt)-(p{Yo,Mo)- f dm<f{Ys,Ms) -aMs + aef E[a{Ys-s', Ms-s')]b{ds') 

Jo 1 Jo 


ds 


- [ dyip{Y,,Ms)ds- [ a{Ys,Ms){<p{Q,Ms) - p}{Ys,Ms))ds, 
JO Jo 


is a local martingale. 
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CRISTOBAL QUININAO 


Let US recall now the Jacob-Shiryaev [71 Theorem 11.2.42 page 86] 

Theorem A.l. There is equivalence between: 

• (YfjMt) is a semimartingale, and it admits the characteristics i.e., {Yt,Mt) writes 

= (ro,Mo)+7W" + B, 

where is the continuous local martingale of the canonical decomposition, B is predictable and 
V is a predictable random measure on IR_|_ x namely the compensator of the random measure 
associated to the jumps of X. 

• For each bounded function g) € the process 


ifiiYu Mt) - V2(Fo, Mo) - [ dyipiYs-,M,_) dBy - [ dmTiYs-,Ms-) dB'T 

JO Jo 

- f {ifiiYt- +y,Mt- +x) - ip(Yt-,Mt-) -ydyip(Yt-,Mt-) - mdmPiYt-, Mt-)} v{ds,dy,dm) 
Jo 


is a local martingale. 

Then, in our case of study, by choosing the characteristics 


= 


[ \l + YsaiYs-,M,-)]ds, Sr = [ \-aMs+ae [ E[aiYs-s', Ms-s')]b{ds') 
Jo ^ Jo ^ Jo 


ds, 


and by 

v{ds,dy,dm) = a{Ys_, Ms_) ds S_Y„_{dy)So(dm), 
we get that {Yt,Mt) is a semi martingale. 


As for the second important Jacob-Shiryaev [71 Theorem III.2.26 page 157] cited in the main text, let 
us now rewrite their general result to our study case. Consider the stochastic differential equation 


(A.2) 


(To, Mo) = i^y,U) 

d{Yt,Mt) = f}{t,Yt,Mt)dt + 5{t,Yt-,Mt-,z){N{du,dt) - q{du,dt)), 


where A/" is a standard Poisson process with intensity measure q{du, dt) = du dt. 

Theorem A.2. Let rj be a suitable initial condition (i.e., a probability on R+), and (3,5 be 

J/3 = (/3^,/3^), a Borel function: R._|_ x —>■ 

= (J^,(5^), a Borel function: R+ x R+ x R^ —>■ R^. 

The set of all solutions to (|A.2ll with initial condition rj is the set of all solutions to a martingale 
problem on the canonical space where the characteristics {B,0,v) are given by 

Bl{w) = f fi'^{s,Ys{w), Ms{w)) ds, v{w,dt x dy dm) = dt x Kt{Yt{w),Mt{w),dy,dm), 


with 


nOC 

Kt{y,m,A) = / l[A\{o}}iSit,u,y,m)) du. 
Jo 


We notice that {Yt,Mt) indeed solves the Martingale problem given by (lA.ll) . therefore it is a solution 
to the equation (IA.2I) for some standard Poisson process, and therefore Lemma |T3] is proved. 
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